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Abstract 

We present model independent high resolution amplitude analyses of CERN data on ir~p —> ii~Tt + n at 
17.2 GeV/c for dipion masses 580-1080 MeV at small momentum transfers on polarized target (Analysis 
I) and unpolarized target (Analysis II). The results for S- and P-wave transversity amplitudes from the 
Analysis I are used in a model independent determination of helicity amplitudes. All three analyses provide 
similar evidence for p°(770) — /o(980) mixing in S'-and P-wave amplitudes. In Analyses I and II the S-w&ve 
amplitudes |<Sd| 2 peak at j o°(770) mass while the P-wave amplitudes \Ld\ 2 dip at /o(980) mass. In both 
analyses the S'-wave amplitudes S T are nearly in phase with longitudinal amplitudes L T and nearly 180° 
out of phase with transverse amplitudes U T indicating resonant phases of p°(770) also in the S- wave for 
both target nucleon transversities r = u,d. The S-wave single flip helicity amplitudes |5i| 2 and the relative 
phases $(£1) — show the same pattern. While the non-flip helicity amplitudes |£o| 2 are small below 

960 MeV, there is a sudden rise above this mass that may reflect the presence of /o(980) resonance. We 
show that the apparent dependence of p°(770) width on helicity in amplitudes |P r | 2 and |iV r | 2 results from 
the interference of transversity amplitudes with definite dipion helicities A = ± 1 and not from violation of 
rotational/Lorentz symmetry. This suggests that the production mechanism and the mechanism responsible 
for p°(770) — /o(980) mixing have independent dynamical origins. 
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I. INTRODUCTION. 



Following the discovery in 1961 of p meson in irN — > ttttN reactions [l[, the measurements of 
forward-backward asymmetry in ir~p — > ir~iT + n suggested the existence of a rho-like resonance in 
the S-wave, later referred to as <r(750) scalar meson [2t6J. It was expected that cr(750) would show 
up prominently in ir~p — > 7r°7r°n production where /?°(770) does not contribute. The measurements 
of this reaction at CERN in 1972 found no evidence for a rho-like <r(750) [7|. Furthermore, in 1973, 
Pennington and Protopopescu used analyticity and unitarity constraints on partial wave amplitudes 
in 7T7T — > tttt scattering (Roy equations) to show that a narrow cr(750) resonance cannot contribute 
to tttt scattering [8]. From these facts it was concluded that <t(750) does not exist in n~p — > ir~-ir + n 
and in 1974 Particle Data Group dropped this state from its listings. 

In 1978 Lutz and Rybicki showed that almost complete amplitude analysis of reactions ttN — > 
ir + Tr~N is possible from measurements in a single experiment on a transversely polarized target [jj. 
CERN-Munich-Cracow group measured ir~p — > ir~ir + n on transversely polarized proton target 
in a high statistics experiment at 17.2 GeV/c at low momentum transfers —t = 0.005 — 0.20 



(GeV/c) 2 [lO-UJ] and at high momentum transfers -t = 0.20 - 1.00 (GeV/c) 2 [13]. Saclay group 



measured tt + u — > Tr + ir~p on transversely polarized deuteron target at 5.98 and 11.85 GeV/c at 
larger momentum transfers —t = 0.20 — 0.40 (GeV/c) 2 [l5|, Qjl]. Measurements of n~p — > 7r~7r + n on 
transversely polarized proton target at 1.78 GeV/c at low momentum transfers —t = 0.005 — 0.20 



(GeV/c) 2 were made at ITEP [17]. 



The CERN measurements of ir p — > ir ir + n and ir + n — > ir + ir p on polarized targets reopened 
the question of the existence of <r(750) scalar meson. Evidence for a narrow <r(750) was found 



in amplitude a naly ses of 7r _ p — > ir~n + n at 17.2 GeV/c [ill, LL3J, LL8|, |20|] and in ir + n — > ir + ir~p at 
5.98 and 11.85 19,|2(j. Clear evidence for <r(750) emerged from later and more precise amplitude 
analyses of both reactions 2114231]. Another evidence for <r(750) comes from the amplitude analysis 
of the ITEP data at 1.78 GeV/c [lj. 

In 2001, E852 Collaboration at BNL reported high statistics measurements n~p — > 7r°7r°n at 
18.3 GeV/c [24]. The data revealed large differences in the S'-wave intensities in 7r~7r + and 7r°7r° 
production. There was no evidence for a rho-like <r(750) resonance in the 7r°7r° S'-wave. The BNL 
data presented anew the puzzle of the a (750) resonance. The CERN and BNL data are both high 
quality data that cannot be used to exclude one another. We must accept that they are both correct 
and that the apparent contradictions between them are telling us something new and important. 
The puzzle of <r(750) resonance has become a unique opportunity to learn new physics. 

The first hints on the solution of the puzzle emerged from the amplitude analysis of CERN 
measurements of ir~p — > ir~ir + n in Ref. 23]. This study extended the dipion mass range of 600 



900 MeV of previous analyses [2l|, |22j to 580 - 1080 MeV to include the Jo (980) resonance. 



The S-wave transversity ampltudes S T with target nucleon trasversity r = u(up) , d(down) were 
fitted using a parametrization with Breit-Wigner amplitudes for both resonances and complex 
backgoudnds. The best fit to all moduli of S-wave transversity amplitudes gives m a = 778 ± 16 
MeV and Y a = 142 ±33 MeV. These values are very close to resonance parameters of p°(770) with 
m p = 775.8 ± 0.5 MeV and T p = 150.3 ±1.6 MeV. Both solutions for the P-wave amplitude with 
zero dipion helicity \Ld\ 2 showed an unexpected dip at ~ 980 MeV, the mass of /o(980) resonance. 

In 2001, Gale and collaborators investigated p°(770) — ao(980) mixing in lepton pair production 
at RHIC arising from the violation of Lorentz symmetry by vacuum but conserved by the interaction 
Lagrangian, and found large measurable effects 25j. Their work suggested that the dip in |Lrf| 2 
arises from / o°(770) — /o(980) mixing in the P-wave. Similar mixing in the S-wave then allows 
to identify <r(750) with / o°(770) which immediately explains why no <r(750) was ever found in 
ir~p — > 7r°7r°n production. The p° (770)- f (980) mixing implies an apparent violation of rotational 
and Lorentz symmetry which appeared to be supported by the observation of large differences in 
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p°(770) widths in the moduli of P-wave amplitues \U T \ 2 and \N T \ 2 with transverse dipion helicities. 

In Ref. [23)] we used the results of the fits to the transversity amplitudes to calculate S-wave 
helicity non-flip and flip amplitudes So and Si assuming the residues of <r(750) were real and 
positive in both transversity amplitudes and there was no additional relative phase apart from the 
fits. The results suggested that <r(750) dominates So while it is suppressed in Si, and therefore 
that it dominates in tttt — > a\TT and is suppressed in tttt — > tttt scattering. However, the results for 
Si disagreed with established tttt scattering amplitudes. In over 100 later fits with different input 
values of complex residues of <r(750) and other fitting parameters we obtained mostly convergent 
fits. They all had the same x 2 and values of resonance parameters of c(750) and /o(980) identical 
to the values from our original fits. However, the corresponding aplitudes So and Si were widely 
different with many of them unphysical solutions. We must conclude that the presence of a rho-like 
resonance in S-wave transversity amplitudes is not related to tttt — > a\TT or tttt — > tttt scattering. It 
stands on its own and lends support to the hypothesis of /9°(770) — /o(980) mixing. 

In this work we present our results of a high resolution amplitude analysis of CERN data on 
polarized target (Analysis I) and of a new model independent amplitude analysis of CERN data 
on unpolarized target (Analysis II). We trace the presence of /3°(770) in the S-wave transversity 
amplitudes in both Analyses to the data component 01 + 02 in the equation |S| 2 = a\ + 02 — 3|L| 2 . 
The data a\ +02 have a large p°(770) peak that is not balanced out by p°(770) peak in the P-wave 
amplitude \L\ 2 . The relative phases between S-wave amplitudes and P-wave amplitudes L T and U T 
are near 0° and —180°, respectively in both Analyses. Since the P-wave amplitudes have p°(770) 
Breit-Wigner phase, the S-wave amplitudes must have nearly the same resonant phase. 

Next we outline a model independent method to determine the S- and P-wave helicity ampli- 
tudes from the knowledge of transversity amplitudes and present new results on helicity amplitudes 
S n and L n , n = 0, 1. The four solutions for the single flip amplitude |Si| 2 show a clear p°(770) 
peak and the phase of Si is again nearly equal to the phase of P-wave amplitude L\. There is a 
remarkable similarity and concordance between the Analyses I, II and the helicity amplitudes. 

Finally we show that the large differences in the widths of \U T \ 2 and \N T \ 2 amplitudes arise 
from the interference of amplitudes of opposite dipion helicities A = ±1, and that the width of 
/9°(770) resonance does not depend on its helicity in both Analyses I and II. This suggests that the 
production process conserves rotational/Lorentz symmetry. This means that the mechanism for the 
p°(779) — /o(980) mixing and the production mechanism have different dynamical origins. While 
there is no reason to assume that the production mechanism is governed by anything other than the 
S-matrix, the p°(770) — /o(980) mixing must originate in a new spin mixing interaction. While the 
production mechanism conserves Lorentz symmetry, the spin mixing mechanism spontaneously 
breaks this symmetry in the measured amplitudes that combine both mechanisms. We address 
these issues in two related papers [2a, [27 ] . 

The paper is organized as follows. In Section II. we define the S-and P-wave subsystem and the 
reduced transversity amplitudes used in all our analyses. In Sections III. and IV. we present the 
analytical solutions for transversity amplitudes with polarized and unpolarized target data which 
define the Analyses I and II, respectively. Evidence for the p°(770) — /o(980) mixing from Analyses 

I and II is presented in Section V. In Section VI. we outline our method of the model independent 
determination of S-and P-wave helicity amplitudes and present evidence for p°(770) — /o(980) 
mixing from amplitudes |Si| 2 and phases &(Li) — <£(Si). In Section VII. we use Analyses I and 

II to demonstrate the helicity independence of the p°(770) width and thus Lorentz symmetry of 
the production mechanism. The paper closes with our conclusions and outlook in Section VIII. In 
the Appendix we show that the measurements of Imp® and Imp® determine uniquely the signs of 
relative phases in the Analyses i and II. 
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II. S- AND P-WAVE SUBSYSTEM OF REDUCED FINAL STATE DENSITY MATRIX. 



The final state density matrix in ir p — > it ir + n has a general form 0, 26] 

p f (0<j>,P) = ~(I°(9<t>,P)a° + T(9<f>,P)a) (2.1) 

where 9, <f> describe the direction of n~ in the dipion centre-of-mass system and P and Q(9(p, P) are 
polarization vectors of the target and recoil nucleon, respectively. The most feasible experiments 
are measurements on unpolarized or polarized targets leaving the recoil nucleon polarization vector 
unobserved. Such measurements provide information on the reduced final state density matrix 
I°(00,P) given by 0,0] 



I°(6<f>, P) = Tr( Pf {e^ ?)) = + PJ^m + P y I y m + P*I°(00) (2.2) 

where the polarization components of I°(9(j),P) are expressed in terms of density matrix elements 

»2 J max 

J<J' A>0 A' 

for k=u (unpolarized target) and k=y (transversely polarized taget normal to the scattering plane), 
and 

j2 J max 

Ilm = iL^ Y^Y,bj>uirn&)$i*wimY$\m ( 2 - 4 ) 

J< J' A>0 A' 

for k=x (transversely polarized target in the scattering plane) and k=z (longitudinally polarized 
target). In (2.3) and (2.4) £ = 1 and £ A = 2 for A > and the factor £jJ' = 1 for J = J' 
and = 2 for J < J' . In a given region of dipion mass m and momentum transfer t only 
amplitudes with J < J ma x contribute and all sums in (2.3) and (2.4) are finite. Prom the measured 
intensity I°(9(j), P) the density matrix elements are determined in small (m, t) bins using maximum 
likelihood method 0, 0, Q . 



Amplitude analyses of measurements on polarized targets are best performed using nucleon 
transversity amplitudes with definite t-channel naturality and N^ T corresponding to unnatural 
and natural exchange, respectively 0, [26|. J and A are dimeson spin and helicity and r = u,d 
is the nucleon transversity for target spin "up" or "down" relative to the scattering plane. In 
ir~p — > 7r~7r + n the amplitudes U^ T exchange tt and a\ quantum numbers in the t-channel while the 
amplitudes -/V^ r exchange 02 quantum numbers. General expressions for density matrix elements 
in terms of transversity amplitudes were tabulated by Lutz and Rybicki and are reproduced in 
the Ref. [26] in our notation. For S- and P-wave transversity amplitudes we shall use a simplified 
notation A T where the amplitudes A = S, L,U, N are defined as 

SV = <r, Lr = Ul T (2.5) 

U T = Ul ;T , N T = Nl T 

Amplitude analysis of the S- and P-wave subsystem of the reduced density matrix determines 
reduced transversity amplitudes defined as follows 
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TABLE I: Density matrix elements Rep^d a/dtdm and Repyd 2 a/dtdm in terms of reduced transversity 
amplitudes. 

p'l'l' Rep a u d 2 a jdidm Rep n y d 2 a / dtdm 

P?s \s? + \s\ 2 \s\ 2 -\s\ 2 

pII \l\ 2 + \l\ 2 \l\ 2 -\l\ 2 

P\\ \{\U? + \U\ 2 ) + i(|Af + \N?) \{\U\ 2 - W\ 2 ) + hi\N\ 2 W\ 2 ) 

pV-i -mi 2 + \m + u\m 2 + \n\ 2 ) ~m\ 2 - p\ 2 ) + m\ 2 - \m 2 ) 



Rep™ 
VZRepH 
V2Rep\° s 


Re(LS* + 
Re(LU* + 
Re(US* + 


LS*) 
LU*) 
US*) 


Re(LS* - 
Re(LU* - 
Re(US* - 


- LS*) 

- LU*) 

- US*) 


P J x{ 


Imp®d 2 a/ 


dtdm 


Imp z d a 


/dtdm 


V2Imp a s \ 
V^ImpH 
Imp]} n 


Re(-SN* - 
Re(—LN* - 
Re(+UN* - 


hNS*) 
^NL*) 
- NU*) 


Im(+SN* 
Im(+LN* 
Im(-UN* 


- NS*) 

- NL*) 
+ NS*) 



S=\S U \, S=\S d \ (2.6) 
L= \L u \expi(<S> Lu -$ Su ), L= \L d \ expi ($ Ld - $ Sd ) 
U = \U u \expi($ Uu - $sj- U = \U d \expi(<S> Ud - <S> Sd ) 

N = \N U \ expi {$ Nu - $ S J , N = \N d \ expi(<S> Nd - $ S J 

where <&a t is the phase of the amplitude A T . The reduced transversity amplitudes are related to 
transversity amplitudes by phase factors 

A u = A exp i$ Su , A d = A exp iui exp i$ Su (2.7) 

for unnatural exchange amplitudes A = S,L,U and 

N u = N exp iu exp i® Su , N d = N exp i<S> Su (2.8) 

for natural exchange amplitude N. In (2.7) and (2.8) <&s u is the arbitrary absolute phase and 
oj = <&s d — &s u is the relative phase between 5-wave amplitudes of opposite transversity. In 
the Table I. we present S- and P-wave density matrix elements expressed in terms of the reduced 
transversity amplitudes. Because of the angular properties of Y^(9(j)), the elements (/^)oo = (Pk)ssi 
(Pfc)oo an< ^ (/°fc)ii' k = u,y are not independent but appear as two independent combinations in the 
measured angular distributions (2.3) 

(pDsp = (p° k )Ts + (pg)So + 2 (P°)n. (pDpp = (p2)oo " (p2)u (2-9) 

At any dipion mass the entire S- and P-wave subsystem can be solved analytically for the 
transversity amplitudes A T up to an overall phase selected to be <J>s u - withount any measure- 
ments of recoil nucleon polarization. The observables Rep^ and Rep® determine the amplitudes 
S, L, U, | TV" | and S, L, U, \N\. The observables Imp® and Imp , determine the phases of amplitudes 
N and N (Appendix). The relative phase uj is determined analytically by the process of conversion 
of transversity amplitudes into helicity amplitudes (Section VII.). The amplitude analysis yields 
two solutions for A u (i),i = 1,2 and two independent solutions for A d (j),j = 1,2 which give rise 
to four different 5-and P-wave final state density matrices Pf(0<f), P;ij). 
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«5 = ^=((p u ) 1 \ + (p y)ll)Z = \L\\U\cOs(<S> LU ) 



III. AMPLITUDE ANALYSIS WITH POLARIZED TARGET DATA ON Rep° u AND Rep° y . 
A. The moduli of transversity amplitudes 

The measured observables p^ and p® organize themselves into two groups involving amplitudes 
of opposite transversity. Using the expressions in Table I. and £ = d 2 a/dtdm, the two groups are 

«i = \{(pI)sp + (p y )sp)Z = \S\ 2 + \L\ 2 + \U\ 2 + |Af (3.1) 
a 2 = ((p u )p P + (p y ) P p)X = 2\L\ 2 -\U\ 2 -\N\ 2 
a 3 = ((p°)^i + (p°)^i)S = |iV| 2 -|f/| 2 

a 4 = ^((p£)S2 + (pS)o2)S = |L||5|co8(* Lfl ) 
1 

a 6 = ^((p£)g + (pS)S)S = |l/||5|co8(^5) 

for reduced transversity amplitudes with transversity t = u, and 

3i = \{(pI)sp ~ (p y )sp)Z = \S\ 2 + \L\ 2 + \U\ 2 + \N\ 2 (3.2) 

«2 = ({pDpp - (p°)pp)S = 2\L\ 2 - \U\ 2 - \N\ 2 

as = ((pS)iii - (p!)iii)S = |iV| 2 - |t/| 2 

«4 = ^(OM - (p° y )l°s)Z = \L\\S\ cos(¥ i5 ) 

«5 = ^((pS)oi " (pS)oi) S = W\ cos(¥ LC7 ) 

a 6 = i=((p°)l°-(p°)^)S = |C7||S|cos(l [/ 5) 

for reduced transversity amplitudes with transversity t = d. The relative phases are defined as in 
(2.6). For dipion masses where S- and P-wave dominate, (Pu)sp and (py)sp are traces 

(p°) 5P = Tr{{pl){{) = 1, (Py)sp = Tr((py)x{) = T (3.3) 

where T is target spin asymmetry. 

The 6 equations (2.11) and (2.12) each involve 7 unknowns for 4 moduli and 3 cosines, and are 
not solvable. The missing equation in each group is supplied not by the data but by phase relations 

$LS ~ $>LU ~ &US = {§L U - $S„) - {<$>L U -&U U )- (®U U - $sj = (3.4) 
$LS - $US - $LU = ($L d - $S d ) - {$L d ~ $>U d ) ~ {®U d ~ §S d ) = 

These conditions lead to non-linear relations between the cosines 
cos 2 (<&ls) + cos 2 ($L(y) + cos 2 (<f>us) - 2cos($L5) 

cos 2 (¥ L 5)+cos 2 (¥ L c/) + cos 2 (¥ [ /5)-2cos(l L 5)cos(l Ll /)cos(¥ l/5 ) = 1 
Substituting for the cosines from (3.1) we get 

al|5| 2 + aj\L\ 2 + a\\U\ 2 - |5| 2 |L| 2 |f/| 2 = 2a 4 a 5 a 6 (3.6) 
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From (3.1) we have three equations for moduli 

'I 2 = (ai 
\U\ 2 = \L\ 2 - ^(a 2 + a 3 



\S\ 2 = ( ai + a 2 ) - 3\L\ 2 (3.7) 
1 



2 

|Af = |L| 2 -i(a 2 -a 3 ) 

Substituting (3.7) into (3.6) we get a cubic equation for \L\ 2 = x 

ax 3 + bx 2 + cx + d = (3.8) 

where a = 3 and 

b = -3[-(aj + 02) + ^(a 2 + 03)] (3.9) 
c = -(ai + a 2 )(a 2 + a 3 ) + af + a 2 , - 3a 2 
d = (ai + a 2 )a 2 - ^(a 2 + a 3 )a| - 2a 4 a 5 a 6 

From (3.2) we obtain equations for moduli \S\ 2 , \U\ 2 , \N\ 2 similar to (3.7) and from (3.5) a cubic 
equation for \L\ 2 = x with observables a^, k = 1, 6 replacing a^, k = 1, 6 in the coefficients (3.9). 

The cubic equations for \L\ 2 and |L| 2 are just another form of the phase conditions (3.4). The 
equations can be solved analyticaly. To solve for x we write x = y — yo and require that the cubic 
equation (3.8) transforms to the form 

y 3 + 3Py + 2Q = (3.10) 
This is accomplished with yo = b/3a and 



fl = sign(Q)yjFf, V=^>0 (3.12) 



Next we define quantities 



There are three categories of solutions of cubic equation (3.10) [16|, [30J. Their numerical study 
using a Monte Carlo method shows there are two physical solutions that have an analytical form 

|L(1)| 2 = - yo + i?cos(|) + \/3i2sin(^) (3.13) 
|L(2)| 2 = - yo + R cos (t)- V3i?sin(^) 

where cos(0) = V, R > and |L(1)| 2 > |L(2)| 2 . It follows from (3.7) that the two physical solutions 
for the moduli IS^i)) 2 , |[/(i)| 2 and |iV(i)| 2 , i = 1,2, have a similar analytical form. Similarly, there 
are two independent solutions for the moduli |^4(j)| 2 , j = 1,2 with R > and |£(1)| 2 > |L(2)| 2 . 
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B. The relative phases of transversity amplitudes 



For each set of solutions for moduli |^4(i)|, i = 1,2 we calculate from (3.1) the cosines 

cosfeW) = wim (3 - 14) 

cos($ Lu(i)) = 



L(i)\\U(i)\ 



\U(i)\\S(i)\ 

with similar equations for cos(&ls(j)), cos($ lu (j)) , cos(& us(j))- From the phase condition 

$ L - $ S = ($L " $£/) + ($u ~ $s) (3.15) 

we obtain 3 equations for sines 

C A = sin(<J» L5 ) sin($ L[/ ) = + cos(<%<?) - cos{<5> LS ) cos{$ LU ) (3.16) 
C B = sin($ L5 ) sin($(75) = + cos($li/) - cos($l,s) cos($[/ 5 ) 
Cc = sinks') sm($ LU ) = - cos($ls) + cos($ us ) cos($ LU ) 



The system is solvable provided that sign ( Cc )= sign (C^Cs). With 

sin($ LS ) = eVl-cos 2 ($ L5 ), (3-17) 

where e = ±1, sin(<J>£[/) and sin^f/s) can be calculated from Ca and Cg, respectively. This 
procedure allows us to determine the phases $t,s(i), ®Lu{i) an d 3>[/s(i). With 

sin(¥ L5 ) = e v l-cos 2 (¥ L5 ) (3.18) 

where e = ±1, we obtain for each set of solutions for the moduli |yl(j)|, j = 1,2 the phases 
&Ls(j), ^w(j), ^us(j)- The amplitude analyses assumed positive root for both sines sin^^s) 
and sin(<l>£s). Such procedure is uniqe provided that the phases $>ls and &ls do not change signs. 
A change of sign would manifest itself as a double zero in the phases <5 ls and ls calculated from 
the positive roots for the sines. We found no evidence of a double zero in the phases §ls an d <&ls 
so calculated, and apart from the sign ambiguity, these phases are uniquely determined. 

The phases <&ls an d &ls are closely related to the positivity of the reduced density matrices 
p°{P y ) = Pu + PyPy f° r P ure initial states P y = ±1 by a relation det(p°(P y )) ~ sin 2 (<I>Ls). The 
fact that these phases do not change signs implies that the eigenvalues of the density matrices 
p°(P y = ±1) are all non-zero. 

For negative roots in (3.17) and (3.18) all phases change signs. There are four combinations 
of the signs e and e. The phases with opposite signs correspond to complex conjugate reduced 
transversity amplitudes. For any given solution of moduli \A(j)\,i,j = 1,2 there is a four- 

fold ambiguity in the phases of the reduced transversity amplitudes A = L, U 

L(i,e) = \L(i)\ eMi^Ls(i)), L(j,t) = \L{j)\ exp(ie¥ L5 (j)) (3.19) 

U(i,e) = \U(i)\ exp(ie$ us (i)), U(j,e) = \U(j)\ exp{ie$us{j)) 

In the Appendix we show that this four-fold sign ambiguity can be resolved into a single set of 
phases for each set of moduli |-A(i)|, \A(j)\,i,j = 1,2 in the process of determination of the phases 
of natural exchange amplitudes from the measurements of Imp® and Imp® on polarized targets. 
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IV. AMPLITUDE ANALYSIS WITH UNPOLARIZED TARGET DATA ON Rep° u . 



We now show that a model independent amplitude analysis is possible with unpolarized target 
data Repu- Perhaps more appropriately, the analysis should be called "intensity analysis". First 
we define partial wave intensities I a = \A U \ 2 + |Ai| 2 = \A\ 2 + |^4| 2 > A = S,L,U,N. Using the 
expressions in the Table I we can define a new set of 6 equations 

(4.1) 



a? = 


\{pI)sp^ = 


= ^(|l5| 2 + |/L| 2 + |%| 2 + |/iv| 2 ) 


a% = 




i(2|/ i | 2 -|/ c/ | 2 -|/^| 2 ) 


4 = 


(p£)iiis = 


\{\In\ 2 -\Iu?) 


4 = 


2(Pu)os^ = 


-|/l||/5|cOs(* L5 ) 


4 = 




= \\Il\\Iu\cos{^ lu ) 


4 = 




= \\Iu\\Is\cos{^ us ) 



where for AB = LS, LU, U S the interference terms are 

\I a \\Ib\ cos(* ab ) = Re(AB*) + Re(AB*) (4.2) 

The system of equations (4.1) can be solved if the correlations cos(^ab) satisfy a cosine condition 
similar to (3.5). To see if such is the case we used (4.2) to expresse the correlations in terms 
of transversity amplitudes. Upon substitution into the cosine condition we obtained a non-linear 
constraint connecting moduli and cosines of relative phases of all transversity amplitudes S, L, U 
and S, L, U. This condition was tested numerically for all four combinations of the amplitudes 
A(i),A(j),i,j = 1, 2 by every single Monte Carlo sampling of the data error volume that produced 
a physical solution. In each case the constraint was satisfied exactly as an identity. We concluded 
that the system (4.1) is analytically solvable using the method described in the previous Section.. 

To prove the analytical solvability in a simpler way we set in (3.1) and (3.2) all observables 
p° y = 0. Then (4.1) reads 

a? = l -{pl) SP T, = \S°\ 2 + |L°| 2 + \U°\ 2 + |iV°| 2 (4.3) 
4 = (p° u )ppZ = 2\L \ 2 -\U \ 2 -\N \ 2 

4 = (p°u)\ 1 ^ = \n°\ 2 -\u \ 2 

4 = \( P ° u )h^ = \L \\S \ cos(d>0 s ) 
4 = ^(p°u) 1 o\^ = \L°\\U°\ co S (^ LU ) 
4 = ^(p°)l°S = |C/°||5 |cos« 5 ) 

where the amplitudes A = A(p^, Py = 0) are analytic continuation of the amplitudes A = A(p^, Py) 
to p°y = 0. We obtain a similar set of equations from (3.2) with a%,k = 1,6 for amplitudes 

A° = A(pu,Py = 0) which are analytical continuation of the amplitudes A(p^,p y ) to p^ = 0. Since 
a£ = a£, k = 1, 6 we must have 

A = A° (4.4) 
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The equations (4.1) and (4.3) coincide when 

l A = \ A \2 + \A\ 2 = 2\A°\ 2 = \A°\ 2 + \A°\ 2 = I° A (4.5) 
cos(^ab) = cos($%) 

Since the amplitudes A and ^4° are complex valued functions, they satisfy the phase conditions 
(3.4) and cosine conditions (3.5). This means that the equations (4.1) are analyticaly solvable. 

Helicity nonflip and single flip amplitudes Ao and A\ are related to transversity amplitudes by 
the relations 0, 26] 

A = -^=(A u + A d ) = ^=(A + Aexp(iu))eM^(S u )) (4.6) 

—i 1 — 

Ai = -j={A u -A d ) = -^(A-Aexp(iu}))exp(i^(S u )) 

where to = &(Sd) — &(S U ). Then at p® = we have 

\A° \ 2 = |A°| 2 + |A°| 2 co S (o;) (4.7) 
|A?| 2 = |^°| 2 -|A°| 2 cos(a;) 

In Section VI. we shall find cos(w) = -1. Thus \A^\ 2 = and \A^\ 2 = 2\A°\ 2 . We thus find another 
useful relation for physical intensities in addition to (4.5) 

I A = \A\ 2 + \A\ 2 = \A° 1 \ 2 = I A (4.8) 

There are two solutions for the single flip helicity amplitude A\ extrapolated to p® = which 
correspond to the two analytical solutions of the equations (4.1) satisfying the conditions (4.5) 

|^(1)| 2 = I A (l,l) = \A(l)\ 2 + \A{l)\ 2 (4.9) 
|A?(2)| 2 = / A (2,2) = |yl(2)| 2 + |^(2)| 2 

The off-diagonal intensities I A (1,2) and 1^(2,1) are solutions of the equations (4.1) that are not 
related to the extrapolation of the amplitudes to p® = and thus do not meet the conditions (4.5). 
The reason for this is the fact that I A (1, 2) ^ I A {2, 1) while I A (1, 2) = I A (2, 1). 

The amplitude analysis determines also the relative phases between single flip amplitudes A®. 
To determine <£(-£?) — ^(S^) we find from (4.6) the interference term 

LxSl = ~(LS* + L5* - I-S* exp(iw) + LS* exp(-iu)) (4.10) 

where S,S are real, L = \L\ exp(z<l>sx) and L = \L\ exp(i§sL)- At p® = we have 5° = S° , 
L° = 1° and <S>° LS = $° LS . Then (4.10) reads 

L?S?* = |L°||5°|exp(i*£ s )(l-cosw) (4.11) 

With cos a; = — 1 we find for both Solutions 1, 2 

$(L?) - $(5?) = $(L°) - $(S?) (4.12) 

where r = u, d. Note how our results on the single flip amplitudes A\ depend on our finding in 
Section VI. that the relative phases $(S d ) - $(S U ) = ±180°. 
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V. EVIDENCE FOR p°(770) - / (980) MIXING: TRANSVERSITY AMPLITUDES. 
A. Data analysis 



We have performed high resolution amplitude analyses of CERN data on ir p — > tt 7r + n at 
17.2 Gev/c at 0.005 < \t\ < 0.20(GeV/c) 2 and in dipion mass range 560 < m < 1080 MeV with 
polarized target data Pu, Py 12, 13] (Analysis I) as well as with only unpolarized target data p° [i^ ] 
(Analysis II). The two analyses were performed using the same code but in the Analysis II we set 
Py = 0. The analysis with polarized target data used 1 million of Monte Carlo samplings of error 
volumes of measured density matrix elements. Not all values of density matrix elements correspond 
to physical amplitudes. The procedure produces frequency distribution of physical solutions for 
the moduli and phases and for various observables calculated from the solutions for amplitudes. 
The average values of the distributions for the moduli add up exactly to X = d 2 a/dtdm and those 
for the phases satisfy cosine conditions (3.5). These averaged values can thus be interpreted as 
the measured amplitudes. The range of distributions determines the asymmetric error bars on the 
amplitudes and calculated observables. Analysis I using 5 million Monte Carlo samplings produced 
virtually identical averaged values and slightly larger errors, indicating the analysis is stable. The 
analysis with unpolarized target data used 200 000 Monte Carlo samplings. 

The results from both analyses are shown in Figures 1-7 and 11-13. The Figures show data 
points with errors from both analyses. 
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FIG. 1: Data components a\ + and a® + a° from Analyses I and II, respectively. 
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FIG. 2: Moduli of P-wave amplitudes \L T \ 2 and |L°| 2 from Analyses I and II (with line). 




FIG. 3: Moduli of S*-wave amplitudes ISVI 2 and |5°| 2 from Analyses I and II (with line). 
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B. Evidence from the moduli 

To understand the origin of the p°(770) — /o(980) mixing in the 5- wave amplitudes we recall 
from (3.7) the relation 

\ S \ 2 = ( ai + a 2 ) - 3\L\ 2 (5.1) 

in the Analysis I and a similar equation for the amplitudes IS" ) 2 in terms of a® + a!] in the Analysis 
II. We show that the presence of /d°(770) in the S'-wave amplitudes in both analyses arises from 
the dominant presence of this resonance in the data components a\ + a 2 and a\ + a[>. The data 
components were calculated from those Monte Carlo samplings of density matrix elements in the 
error volume of the data for which physical solutions for the amplitudes were found. The results 
for ai + d2 and a? + a[> are shown in Figure 1. The Figure shows pronounced p°(770) peak for 
target spin "down" component and a suppression of /9°(770) in the target spin "up" component. 
There are pronounced dips at 970 MeV and 1030 MeV in target spin "down" and "up" components, 
respectively, corresponding to the _fo(980) resonance. In contrast, the components a^+aS] are e< l ua l 
for both target transversities, peak at p°(770) mass and dip at 980 MeV. 



Solution (1,1) Solution (2,2) 
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FIG. 4: Intensities I a and I A for S- and P-wave amplitudes from Analyses I and II (with line). 

Figure 2 shows the two solutions for the P-wave amplitudes \L U \ 2 and \Ld\ 2 from the Analysis 
I and the amplitudes |P^| 2 = \L®\ 2 from the Analysis II. Remarkably, in both analyses the two 
solutions of the two cubic equations reproduce closely the general features of the data components 
in the Figure 1. In the Analysis I the p°(770) peak is enhanced in \Ld\ 2 and suppressed in |P U | 2 . 

Figure 3 shows the two solutions for the S'-wave amplitudes \S U \ 2 and \Sd\ 2 from the Analysis I 
and amplitudes |<S^| 2 = |S^| 2 from the Analysis II. The rho-like resonance and /o(980) are clearly 
resolved in both solutions for \Sd\ 2 in Analysis I. This pronounced peak at 770 MeV with a width 
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at half-height of ~ 155 MeV has its origin in the /9°(770) peak of the data component a\ + a>2 which 
peak survives the subtraction of 3|L^| 2 . In the Analysis I the p°(770) peak is enhanced in \Sd\ 2 
and suppressed in \S U \ 2 . The Analysis II shows a clear p°(770) peak in Solution 1 and a subdued 
p°(770) peak in Solution 2 in the amplitude |£^| 2 . 

In both analyses the amplitudes \Sd\ 2 and \S®\ 2 show a shoulder at ~ 930 MeV in Solution 
1 and a secondary peak at the same mass in Solution 2. In Ref. [23| we have shown in fits to 
the moduli of S- wave amplitudes that these structures arise from the interference of p°(770) and 
/o(980) Breit-Wigner amplitudes with a complex backgroud. 

In Figure 4 we present the partial wave intensities Is{hj)^Lihj) from the Analysis I and 
j), j) from the Analysis II for the Solutions = (1,1), (2, 2). The intensities from 

the two analyses are virtually identical, although they come from two different calculations with 
differenet sets and numbers of physical Monte Carlo samplings. This result confirms the validity 
of our analysis with unpolarized target data. Both analyses show the expected dominant p°(770) 
peak in the Solutions I<?(1, 1) and 1). The crucial result however is the observation that both 
analyses show a clear albeit subdued p°(770) peak in the Solutions 1^(2, 2) and IgC2, 2). This means 
that an exact amplitude analysis of unpolarized target data provides evidence for the presence of 
p°(770) in both single flip amplitudes |5?(1)| 2 and \Sf(2)\ 2 . 

The observation of p°(770) — /o(980) mixing in the S-wave amplitudes brings up the question of 
/ o°(770) — /o(980) mixing in the P-wave amplitudes. Figure 2 shows a clear dip at 970 MeV in the 
amplitude \Ld\ 2 followed by a sharp rise. The dip occurs at the mass of scalar resonance /o(980) 
and is more pronounced in the Solution 2. Figure 11 shows unexpected resonant structure (a bump) 
at the /o(980) mass in the Solution 1 of the amplitude \Ud\ 2 ■ These findings are a strong indication 
for p (770) — /o(980) mixing in the P-wave amplitudes. Since /o(980) resonance is produced by 
unnatural exchange we expect /o(980) to be suppressed in natural exchange amplitudes |iV T | 2 as 
evidenced by the data in Figure 12. 

C. Evidence from the relative phases 

Figures 5, 6 and 7 show the mass dependence of the phases <&ab and &ab, AB = LS,US, LU 
from the Analysis I. As the figures show, these phases are continous and nearly constant functions. 
From this fact we can understand why the phases &lSi&ls cannot change sign. From (3.16) we 
see that a change of signs of &ls or results in the change of sign of all phases of the same 
transversity. This would result in large unphysical discontinuities in the phases &uS>&US an d 

Figures 5, 6 and 7 also show the mass dependence of the phases $\ B , $ab> AB = LS,US, LU 
from the Analysis II. There is a remarkable agreement between the Analyses I and II except for 
Solution 1 above 900 MeV in phases $° LS , ® LS , <3?^ 5 , <J>^ 5 . While the phases §ab and ®ab are 
similar, the phases and &ab are equal. Due to the similarity of phases §ls and &ls we 
may conclude that not only $(L?) - $(5?) = $(L°) - $(S°),r = u,d but also $(L?) - $(5?) « 
$(L r ) - $(S T ) for both Solutions 1 and 2. 

Independent evidence for the presence of /9°(770) in the S-wave amplitudes comes from the 
relative phases &ls = &L U — &s u and &ls = &L d — &s d from both analyses. For target spin "up" 
the phase &ls is constant and small, reflecting the similarity of \S U \ 2 and \L U \ 2 in both Solutions 
1 and 2. For target spin "down" the phase $^5 is near zero in the Solution 1 and small and slowly 
varying in Solution 2 below 1000 MeV. These results indicate that the amplitudes S T and L T are 
nearly in phase. Since the amplitudes L T resonate at p (770), so must the amplitudes S T . 

This conclusion is in agreement with the phases <&us = $u u — &s u , &US = &u d — &s d from 
both analyses. These phases are approximately constant below 1000 MeV and indicate that the 
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FIG. 5: Relative phases $ls = $l u — 3>s„ an d ®ls = $L d — $s d from Analyses I and II (with line). 
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FIG. 7: Relative phases <§>lu = $l„ — 3>t/„ and <& LU = § Ld - $ Ud from Analyses I and II (with line). 

amplitudes S T and U T are nearly 180° out of phase. Since amplitudes U T resonate at p°(770), so 
must the amplitudes S T . The Figure 7 shows that the relative phases of the amplitudes L T and U T 
are nearly constant and nearly 180° out of phase, as expected from these two P-wave amplitudes. 

VI. EVIDENCE FOR p°(770) / (980) MIXING: HELICITY AMPLITUDES. 
A. Helicity amplitudes 

Amplitude analysis of the complete S- and P-wave subsystem of the reduced density matrix 
yields four sets of solutions for reduced transversity amplitudes A(i),A(j),i,j = 1, 2. For unnatural 
exchange amplitudes A = S, L,U the reduced transversity amplitudes are related to transversity 
amplitudes by phase factors 

A u (i) = A(i)e W i<S> Su {i) (6.1) 
Ad(j) = A(j)ex.piuJi j expi<S>s u ('i) 

where $5„(i) is the arbitrary absolute phase and 

Wv = *S d 0')-*S u (i) (6-2) 

is the relative phase between S'-wave amplitudes of opposite transversity. Recall from (2.6) that 
S(i),S(j) are real and positive amplitudes, the phases of L(i) and L(j) are $>L,s(i) an d &Ls(j), 
and the phases of U(i) and U(j) are ^us(i) and &us(j), respectively. 
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Helicity amplitudes A J Xx Qu with definite i-channel naturality were defined and related to 
transversity amplitudes of definite i-channel naturality in Ref. 0, [2(| for any dipion spin J and 
helicity A. Due to the P-parity conservation only helicity nonflip and helicity flip amplitudes 
A J X q = A J X+ 0+ and A x x = A J X+ _ are independent. Here n = 0, 1 is nucleon helicity flip n = \\— v\ ■ 
The S- and P-wave helicity amplitudes A n are related to the transversity amplitudes A T by rela- 
tions im 

(-i) n 

A n = K —jL-{A u + {-l) n A d ) (6.3) 

In terms of reduced transversity amplitudes we can write for the unnatural exchange amplitudes 

(—j) n — 
A n (ij) = ^=-(A(i) + (-l)M(j) expCiwy)) exp(i<&^(i)) (6.4) 

Each set of solutions comes with a fourfold sign ambiguity in the phases of the amplitudes 

(A(i), A(j)) e e where ee = ++, H — , — h, are the signs of relative phases &Ls(i) and $>lsO)- 

From (3.16) we see that the change of sign of $Ls(i) (&lsU)) results in the change of sign of the 
phase of amplitude U{i) (U(j)), or complex conjugation of amplitudes A(i) (A(j)). We can thus 
write the four sets of phases for each i,j 

(A(i),A(j)) ++ = (A(i),A(j)) (6.5) 
(A(0,A0'))+- = (A^Aij)*) 

= (A(i)\A(j)) = ((A( l ),AU)) + -T 

(A(i),A(j))^ = (A(iy,A( j y) = ((A( i ),A( j )) ++ y 

As shown in the Appendix, the measurements with planar target polarization can select uniquely 
one of the four solutions for the phases. Since the transversity amplitudes A u (i) and A^(j) are also 
complex conjugated, we get for helicity amplitudes 

A n {i3)- + = {-) n {A n {ij) + _y (6.6) 
Aniij)- = (-) n (A n (ij) ++ y 

Since we do not get new solutions with the phase sets and — + numerical calculations were 

done only for the sets with phases +- 1- and H — . 



B. Bilinear terms of the helicity amplitudes 

We now look at the bilinear terms A n A* = |A n | 2 , A = S,L,U and A n P*,AP = LS,US,UL. 
For the sake of brevity we shall omit in the following the indices ij and ++, H — . Using (6.4) or 
(6.6) we obtain 

|A„| 2 = ^(|A| 2 + |A| 2 + (-ir2X A cos(a;) + (-ir2y A sinH) (6.7) 
= \{lA + (-l) n 2X4Cos(w) + (-l) n 2y A sinH) 

where Xa = Re(AA*), Ya = Im(AA ) and I a = |A| 2 + |A| 2 is partial wave intensity. Note that 
Ys = Im(SS ) = as both S and S are real. For the bilinear terms A n B* we obtain 

A n B* n = ^(aB* + AB* + (-l) n (APV^ + AB*e +iu) )^j (6.8) 
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The real part reads 



Re(A n B* n ) = ^(Re(AB*) + Re{AB*) 



(6.9) 



+(-l) n ((Re( y AB*) + Re(AB*)) cos u + (Im(AB*) - Im(AB*)) sin w)) 
The imaginary part reads 

Im{A n B* n ) = - (lm{AB*) + Im(AB*) (6.10) 

+ (-l) n ((/m(AB*) + /m(AB*))cosu; - (i?e(AB*) - i?e(AB*)) sin w)) 

It is apparent from (6.4) and (6.6) that the knowledge of u> allows to determine the helicity 
amplitudes up to an absolute phase for each solution set i, j and phase set ++ and H — . As we 
show in the next Section, the phase ui can be determined analytically from the consistency condition 

\A n \ 2 \B n \ 2 = (Re(A n B*)) 2 + (Im(A n B* n )f (6.11) 

where the terms are given by (6. 7), (6. 9) and (6.10). 



C. Analytical solutions of the relative phase cu 

In order to make use of the consistency condition (6.11) to determine oj, we first rewrite the 
terms (6.8), (6.10) and (6.11) in a simplified form. To this end we define 

X(AB) = Re(AB*) + Re(AB*) (6.12) 

Y(AB) = Im(AB*) + Im(AB*) 
X(AB) + = Re(AB*) + Re(AB*) 
Y(AB) + = Im(AB~*) + Im(AB*) 
X(AB)_ = Re(AB*) - Re(AB*) 
Y(AB)_ = Im(AB~*) - Im(AB*) 

Then 

Re(A n B*) = ^(X(AB) + (-l) n (X(AB)+ cos w + Y(AB)_ sinw)) (6.13) 
Im{A n B* n ) = ^(Y(AB) + (-l) n (Y(AB) + cosuj - X(AB)-smoj)^ 
\A n \ 2 = ^X(AA) + ( y -l) n (X(AA) + cosuj + Y(AA)-smuj)^ 
Note that X(AA) = I A , X(AA)+ = 2X A and Y(AA)_ = 2Y A . We can write 

Re(A n B* n ) = l -(x{AB) + (-l) n XG{AB)) (6.14) 
Im(A n B* n ) = l -(Y(AB) + (-l) n YG(AB)) 
\A n \ 2 = l -(x{AA) + {-l) n XG{AA)) 
l -(x(BB) + (-lTXG(BB)) 



\B n \ 2 
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where 

XG(AB) = X (AB) + cos u + Y(AS)- sin u (6.15) 
YG(AB) = Y \AB) + cos u - X(AB)_ sin u 

Next we require that 

\A n \ 2 \B n \ 2 = {Re{A n B* n )f + {Im{A n Bl)f (6.16) 

The l.h.s. of (6.16) reads 

X(AA)X(BB) + XG(AA)XG(BB) (6.17) 

+{-l) n X(AA)XG(AA)X + (-l) n X(BB)XG(AA) 
The r.h.s. of (6.16) reads 

X{ABf + XG(AB) 2 + Y(AB) 2 + yG(AB) 2 (6.18) 

+(-l) n 2X (AB)XG(AB) + (-l) n 2F(AB)yG(AB) 

Subtracting (6.16) with n = 1 from (6.16) with n = and using (6.15) we obtain equation linear 
in cosw and sinw 

sin u (x{AA)Y{BE)- + X(BB)Y(AA)_ - 2X{AB)Y{AB)_ + 2Y(AB)X{AB)^j = (6.19) 

- cos uj (x (A A) X{BB)+ +X(BB)X(AA)+ - 2X{AB)X(AB) + - 2Y(AB)Y(AB) + ^j 
which can be cast in the form 

sin u W 2 = - cos wWi (6.20) 

Using sin 2 oj + cos 2 oj = 1 we find 

±W * f(t ou 

coslo = (6.21) 

TU 

SmW = ~lu~ 

where W = ^/Wf+W^. 

Using (6.7), (6.9) and (6.10) it is straightforward to verify that the consistency condition (6.16) 
reduces to identity in the following two cases 

cosw = 0, sinw = ±1 (6.22) 

sinw = 0, cosoj = ±1 (6.23) 

The first case leads to two solutions oj = +tt/2 and oj = —tt/2. The second case leads to another 
two solutions oj = and uj = ir. 
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FIG. 8: Moduli \So(ij)\ 2 and \Si(ij)\ 2 for i,j = 1,2. The solutions for the phase signs ++ and +-arc equal. 



D. Solutions for S- and P-wave helicity amplitudes 

1. Numerical calculations and their checks 



The Monte Carlo amplitude analysis of reduced transversity amplitudes and related observables 
was carried out by a computer code A. The mean values of moduli and phases of the reduced 
transversity amplitudes satisfy strict normalization and phase conditions, respectively, and thus 
represent the true measured amplitudes at the t-bin average of 0.067 (GeV/c) 2 . These mean values 
were used as an input in an exploratory computer code B to calculate to using (6.21) and helicity 
amplitudes A n ,n = 0, 1, A = S, L, U using uj from (6.21), (6.22) and (6.23). Errors on uj from 
(6.21) and on helicity amplitudes were not calculated due to high non-linearity of the equations. 

Equations (6.21) were used with amplitudes A = L, B = S to calculate uj(ij) for each set ++ 
and H — of signs of phases. Using these values of uj(ij) as well as the values corresponding to the 
solutions (6.22) and (6.23), moduli |A n | 2 , \B n \ 2 and interference terms Re(A n B^) , Im(A n B*) were 
calculated to determine cosines and sines of the relative phases 

<$>{A n B* n ) = <S> An -<S> Bn (6.24) 

Since the calculations of co, moduli and the interference terms are all entirely independent, the 
selfconsistency of the data and calculations was checked using the following three tests on the 
relative phases. The first test are trigonometric identities 

cos 2 <5>(A n B*) + sin 2 $(A n B*) = 1 (6.25) 
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FIG. 9: Moduli \L (ij)\ 2 and \Li(ij)\ 2 for ij = 11,22. The solutions for the phase signs ++ and +- are 
different. 

The second test are phase conditions 

(*L n ~ *sj - (<%„ - 4>sJ + (<&t/ n - $ lJ = (6.26) 

The third test are cosine conditions equivalent to phase conditions 

cos 2 $(L n S; t ) + cos 2 $(?7 n 5;) + cos 2 $(?7 n L;)-2cos $(L„S*) cos &(U n S*) cos = 1 (6.27) 

All conditions are satisfied identically within the single precision calculation used by the code for 
all combinations of solutions i,j = 1,2 for both sets ++ and H — of signs of phases. For instance, 
typical deviation from 1 of the trigonometric identities (6.25) are of order 10 -5 — 10~ 7 . 

To select a physical solution we used two criteria: (1) The single flip amplitudes must dominate 
non-flip amplitudes due to the pion exchange dominance at small t. (2) The amplitudes L\ and 
U\ must show a clear p°(770) peak. 

The computer codes A and B were merged in a computer code C to calulate errors on helicity 
amplitudes using a physical solution for uj found in the exploratory code B. Each Monte Carlo 
selection from the data error volume that yields a physical solution of transversity amplitudes also 
yields the corresponding physical solution of helicity amplitudes, enabling code C to calculate their 
average value and error range in each mass bin in the same way as for the transversity amplitudes. 
The transversity and helicity amplitudes with errors thus represent identical experimental data. 



2. Unphysical solutions for helicity amplitudes 

The values of and uj(ij)-\ calculated from (6.21) show random variations as a function 

of dip ion mass for all solution sets i,j = 1,2 and signs of phases. All amplitudes exhibit the 
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same random behaviour and do not show the required resonant Breit-Wigner behaviour at p°(770) 
resonance. Moreover, the non-flip amplitude Lq has magnitude comparable to or larger than the 
flip amplitude L\. These solutions with cosw / and sin uj / are rejected as unphysical. 

The two values for sino; = ±1 from (6.22) correspond to = ±7r/2 and w^-^ = ±7r/2 

for all solution sets i,j = 1,2. The solutions with sinw = +1 are excluded because they require 
that the non-flip amplitude Lq is larger than the flip amplitude L\. Because cosw = and the 
amplitudes So and S\ do not depend on sinw, both solutions sinw = ±1 give |So| 2 = |Si| 2 = Is/2 
where Is is S-wave intensity. As a result, both values sino; = ±1 are unphysical for any solution 
set i,j = 1, 2 and signs of phases ++ and H — . 

The solutions with cosw = +1 (uj = 0) from (6.23) require that (Sil 2 < \Sq\ 2 and |Li| 2 < |L | 2 
for any solution set i,j = 1,2 and signs of phases. The magnitudes of Si are small compared to 
magnitudes of So, in contradiction with the pion exchange dominance. These solutions are thus 
excluded as unphysical solutions. 



3. Physical solutions for amplitudes S n and L n 

In the solutions with cos a; = — 1 [u = ±180°) the pion exchange dominance is observed in all 
amplitudes and |-Li| 2 and \Ui\ 2 show a well defined p°(770) peak. The resulting amplitudes thus 
represent a unique solution to both helicity and transversity amplitudes, up to signs of the phases 
and to be resolved by the measurements of Im(P x and Imp®. The solution cos a; = — 1 
was used in the computer code C to calulate average values and errors of helicity amplitudes in a 
Monte Carlo analysis of the data. 



23 



Figure 8 shows the moduli of S- wave helicity amplitudes. The solutions for the phase signs 
++ and H — are equal. All solutions for \Si(ij)\ 2 , i, j = 1,2 provide evidence for p°(770) — /o(980) 
mixing in the S- wave helicity flip amplitudes. The p°(770) peak stands out in in Solutions |Si(ll)| 2 
and |Si(21)| 2 which resemble the Solution 1 of \Sd\ 2 - The peak is clearly visible in the Solutions 
|Si(12)| 2 and |Si(22)| 2 which are similar to Solution 2 of \Sd\ 2 - The secondary structure in the 
region of 900-1000 MeV is due to Jo (980). 

Figure 9 shows the moduli of helicity amplitudes L and L\. The Solutions (2, 1) and (1, 2) (not 
shown) are similar to Solutions (1,1) and (2,2), respectively. The helicity flip amplitudes |Li| 2 
dominate the small non- flip amlitudes with the p°(770) peak below 980 MeV. Above 980 MeV 
|Li| 2 decreases rapidly while there is a sudden rise in |Lo| 2 - These structures are associated with 
/o(980) resonance and support the suggestion of p°(779) — /o(980) mixing in the helicity non-flip 
amplitude Lq. The exception is the Solution (2, 2) + + where the rise is observed in |-Li| 2 . 

Figure 10 shows the relative phase ^(^i) — 3>(<Si) which provide an independent evidence for 
p°(779) — /o(980) mixing in the S-wave. Again, the Solutions (2,1) and (1,2) (not shown) are 
similar to Solutions (1,1) and (2,2), respectively. In all solutions the phase is small and nearly 
constant below 960 MeV. Since L\ resonates near p°(770) with a Breit-Wigner phase ^(p ), the 
phase of the amplitude Si must be nearly the same resonant phase. The sudden rise of the phase 
^(Li) — <&(Si) in the /o(980) mass region is due to the presence /o(980) resonance in the amplitude 
Si. 

The non-flip amplitudes |So(ij)| 2 are very small and nearly constant. As a result we find 
\Si(ij)\ 2 ~ Is(i,j) for all i,j = 1,2. Using (4.9) this means that 

|S°| 2 » |Si| 2 (6.28) 

There is also a clear similarity between the phase $(ii) — $(Si) and the phases &ls an d &ls 
which implies that 

- $(S?) « $(Li) - $(5i) (6.29) 

In equations (6.28) and (6.29) the Solutions i = 1, 2 on the l.h.s. correspond to Solutions ij = 11, 22 
on the r.h.s., respectively. The extrapolation of Si to p® = obtained in the analysis of unpolarazied 
target data carries nearly as much information as the exact amplitude from analysis of polarized 
target data. 

E. A note on tttt scattering amplitudes 

The residues of the single flip helicity amplitudes Si and L\ at the pion pole t = p 2 in the 
unphysical region are related to the S-wave and P-wave partial wave amplitudes fs and fp in 
7r~7r + — > 7r~7r + scattering. Analyticity and unitarity of tttt scattering amplitudes (Roy equations) 
exclude a rho-like resonance in the S-wave [§]. Past extrapolations of unpolarized target data to 
the pion pole produced a residue of |S° (2, t = p 2 , m)| 2 with a broad structure in dipion mass below 
900 MeV without any evidence for a rho-like resonance while the residue of \Si(l,t = p 2 ,m)\ 2 
showed the p°(770)-like peak. The question arises whether the extrapolation of our single flip 
amplitudes |Sj*(2, t, m)\ 2 or |Si(12,t,m)| 2 and |Si(22, t, m)\ 2 to the pion pole preserves the p°(770) 
peak or smears it out. The broad shapes of moduli |S?(2)| 2 , |Si(12)| 2 and |Si(22)| 2 are all very 
similar and strongly suggest that an extrapolation of these Solutions to the pion pole will show 
a broad structure consistent with previous extrapolations of unpolarized data. In this regard our 
results agree with the previous tttt scattering data analyses. 

Roy equations only inform us that we should use the extrapolation of |Sj ) (2, t, m)\ 2 and not 
that of |Si(l,t, m)\ 2 , or extrapolations of |Si(12, t, m)\ 2 and |Si(22, t, m)\ 2 and not those of 
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\Si(U,t,m)\ 2 and |Si(21, t, m)| 2 , to calculate tttt scattering amplitudes. Roy equations have noth- 
ing to say about the physical amplitudes and intensities in n~p — > ■n~TT + n scattering process itself. 
Roy equations thus cannot be used to select and reject solutions for physical intensities and ampli- 
tudes in the tr~p — > ir~7r + n scattering, and so deny any physical meaning to the multiple solutions 
for amplitudes which may have - and actually do have - such a meaning apart from and outside of 
the 7T7r scattering theory. 



VII. p°(770) - Jo (980) MIXING AND LORENTZ SYMMETRY. 

It is generally expected that the position and the width of the p°(770) peak as observed in 
the spin averaged cross-section will be faithfully reproduced on the level of spin amplitudes. The 
CERN data on polarized target show that this is not the case. From Figures 2, 11 and 12 we see 
that the p°(770) production is suppressed in all target spin "up" amplitudes while the target spin 
"down" spectra dominate the /3°(770) production. The width at half- height of the peaks in the 
longitudinal spectra \L U \ 2 and | -ZL-rf | 2 is the expected ~ 150 MeV. However, the p°(770) width shows 
different values in different transverse spectra. In the spin "down" spectra the p°(770) width is 
narrower at ~ 120 MeV in \U d \ 2 but wider at ~ 180 MeV in |iV d | 2 . These values are reversed in the 
spin "up" spectra with ~ 180 MeV and ~ 120 MeV in amplitudes \U U \ 2 and |iV u | 2 , respectively. 
Such large variations in the p°(770) width are entirely unexpected and appear anomalous. 

The rotational symmetry of strong interactions requires that the width and mass of a resonance 
do not depend on its helicity. It also prevents the mixing of scalar and vector resonances in the same 
partial wave with a definite spin J. If rotational symmetry was broken in pion production then 
the widths of p°(770) in the P -wave amplitudes L T , U T and N T could differ and p°(770) — /o(980) 
mixing could occur in both S- and P-wave amplitudes. 

But how realistic is the hypothesis that rotational invariance is violated in pion production? 
The breaking of rotational symmetry implies a breaking of Lorentz symmetry. In 2001, Gale 
and collaborators examined vector - scalar mixing of p (770) — ao(980) resonances in dilepton 
production in hot and dense matter. The medium-induced breaking of rotational and Lorentz 
symmetry leads to mixing of different spin states even when the interaction Lagrangian respects 
all required symmetry properties [25| • Their work inspired the idea to test rotational and Lorentz 
symmetry in CERN data on n~p —> Tr~ir + n on polarized target. 

To test the rotational symmetry we need information on transversity amplitudes with def- 
inite dipion helicity A = 0, ±1. The transverse amplitudes U T and N T are a mix of transverse 
amplitudes with dipion helicities A=+l and -1, and thus are not suitable to test the rotational 
symmetry. The required amplitudes with transverse helicities and H~ l are related to the 
amplitudes U T and N T 0, 26] 

H+ 1 = ^=(U T + N T ) (7.1) 

h~ 1 = ±={u u -n u ) 

Their partial wave intensities and polarizations can be calculated from the data on polarized target 

I(H T ) = \H+ 1 \ 2 + \H- X \ 2 = \U T \ 2 + |iV r | 2 (7.2) 
P(H T ) = 2Re{H+ l H- 1 *) = \U T \ 2 - \N T \ 2 

For amplitudes with zero helicity we have longitudinal amplitudes L T . It is convenient to relabel 
them as H® = L T . Their moduli squared shown in Figure 2 are then the longitudinal intensities 

I{H° T ) = \H° T \ 2 (7.3) 
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FIG. 11: Moduli of P-wave amplitudes \U T \ 2 and \U°\ 2 from Analyses I and II (with line). 




FIG. 12: Moduli of P-wave amplitudes |7V T | 2 and \N°\ 2 from Analyses I and II (with line). 
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FIG. 13: Solutions (1,1) and (2,2) for the intensities I{H d ) and I°(H d ) from Analy scs I and II. 
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Figures 13 and 14 show the transverse intensities I{Hd) and I(H U ) for two solutions (1,1) and 
(2,2). The intensities show a clear peak with the same width at half-height of ~ 150 MeV for both 
target spins in both solutions. This transverse p°(770) width is exactly the same as the longitudinal 
width at half-height in the intensities I(H^) and I(H®). This indicates that the p°(770) poles in 
all helicity amplitudes H^, A = 0, ±1 have the same width ~ 150 MeV. 

To explain the large differences in the p°(770) witdth observed in the spectra \Ud\ 2 and \Nd\ 2 
we note from (7.2) 

\U T \ 2 = 1 -{I{H T ) + P{H T )) (7.4) 
|iV r | 2 = 1 -{I{H T ) - P{H T )) 

The calculated polarization P(H c i) is negative and broad in the p°(770) mass region. As the result 
of the opposite signs of P{H d ) in (7.4), the /3°(770) width is narrower in \Ud\ 2 and wider in |-A^| 2 . 
The calculated polarization P{H U ) is positive and broad around the p°(770) mass and therefore 
it has an opposite effect on their widths, as observed. Thus the apparent difference in resonance 
widths in the amplitudes \U T \ 2 and \N T \ 2 is entirely due to the interference of the amplitudes H^ 1 
and R d x . 

Figures 11 and 12 show also the results for amplitudes |C/°| 2 and |iV^| 2 from the Analysis II. It 
is remarkable that in this analysis the resonance widths at half-height in both amplitudes and in 
both Solutions 1,2 are the same at ~ 150 MeV. From (7.4) this implies P°(H T ) = \U°\ 2 - \ N®\ 2 « 0, 
as observed in Figures 11 and 12. The same widths are therefore observed also in both solutions 
(1,1) and (2,2) for the transverse intensities I°(H d ) = \U%\ 2 + \N%\ 2 = I°(H U ) in Figures 13 and 
14. 

We conclude from these tests in the Analyses I and II that the production mechanism con- 
serves rotational/Lorentz symmetry. In contrast, the spin mixing mechanism responsible for 
p°(770) — /o(980) mixing spontaneosly violates this symmetry in the measured amplitudes. While 
the production mechanism is governed by the S'-matrix dynamics of the Standard Model, the spin 
mixing mechanism must be governed by a different dynamics outside of the Standard Model. 

VIII. CONCLUSIONS AND OUTLOOK. 

We have presented experimental evidence for p°(770) — /o(980) mixing in S- and P-wave 
transversity amplitudes in ir~p —> 7r~7r + n from both polarized and unpolarized target data mea- 
sured at CERN. A model independent determination of S"-and P-wave helicity amplitudes confirms 
the evidence from transversity amplitudes and specifies the relative phase of amplitudes with op- 
posite transversities. All three analyses yield mutually consistent results. We show that the 
measurements of Imp® and Imp® resolve the sign ambiguity of relative phases. As a result, there 
are four sets of transversity amplitudes and four sets of helicity amplitudes. 

The principal difference between the Analyses I. and II. is the spin dependence of the moduli 
|>1 T | 2 . In the Analysis I. the spin up moduli \A U \ 2 are suppressed compared to the spin down 
moduli \A(i\ 2 while in the analysis II. these moduli are equal \A U \ 2 = {A^l 2 - The suppression of 
\A U \ 2 in the Analysis I. is a part of a more complex spin dependence of the amplitudes \A T | 2 which 
show oscillations as a function of momentum transfer t in the / o°(770) mass region 

The amplitude Analyses I and II provide the first direct experimental evidence for the inde- 
pendence of the width of a resonance on its helicity expected theoretically from the rotational 
symmetry of strong interactions. Since the interactions in the Standard Model do not violate 
Lorentz symmetry then the observed p°(770) — /o(980) mixing must arise from a different kind 
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271 ] we propose to 



of interaction involved in the pion creation process. In two related papers [ 
identify this new kind of interaction with a non-unitary interaction of the produced 5"-matrix final 
state Pf(S) with a quantum environment E. The four sets of transversity amplitudes and the four 
sets of helicity amplitudes provide two equivalent descriptions of this non-unitary interaction for 
dipion masses below 1080 MeV. 



Determination of the phases of natural exchange amplitudes and the resolution of sign 
ambiguities in the phases $5^ and 

The S- and P-wave density matrix elements Imp® and Imp® involve interferences between 
reduced unnatural and natural exchange amplitudes of opposite transversity shown in Table I.. 
First we shall show that the elements Im{p° x )®\ and Im(p° z )®\ determine the phases a at = &N u — &s d 
and ~a.N = <E>7v d — &s u °f the natural exchange amplitudes N and N, respectively, with a two- fold 
ambiguity for each set of moduli |iV(i)|, \N(j)\,i, j,= 1,2. From the Table I. and (2.6) we have 

ri = >/2Im(p°)£E = Re(-\S\N* + N\S\) = -\S\\N\ cosa N + \N\\S\ cosa N (1) 

si = V2Im{p° z )° s {Z = Im(+\S\N* - N\S\) = -\S\\N\sma N - \N\\S\ sinew 
where £ = d 2 a jdmdt. From (1) we can solve for cos a at and sinojjv 

cosoat = — — (r x + |5||JV|cosajy) (2) 

sinav = — (s\ + 1511 AH sinajv 

|iV||5| V 

Substituting into cos 2 a at + sin 2 ajy = 1 we find 

cos a,N = = — ( A — 2\S\\N\si sinajyl (3) 

2|5||JV|n V 1 " 1 J w 

where 

A=\N\ 2 \S\ 2 - \S\ 2 \N\ 2 - (r 2 + s 2 ) 

Substituting into cos 2 + sin 2 a at = 1 yields a quadratic equation for sin a 7 at with two solutions 

A 



8111 OtN 



2\S\\N\(r 2 + s 2 ) 
where 

4|S| 2 |iV| 2 (r 2 + S 2 ) 



1 ±nVB- 1 (4) 



B 



A 2 

From (4) we can now calculate cosoat 

C ° S ^ = 2|5||iV|fr 2 + S 2 )( ri T SlVWZrt ) (5) 

and from (3) we obtain 

cos aw = - —7=77 9 ~-[r\C ^ s\y/B — 1) (6) 
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where 

ri _ A + 2{rl + s\) 
C ~ A 

Note that for each set of moduli |^4(i)|, |^4(j)|, i,j = 1,2 the two solutions for the phases 
a N(ij, c), a]\f(ij, a) depend i,j and the sign a = ±1 in (4) but do not depend the sign ambi- 
guities e,e of amplitudes A = L, U in (3.19). As a result, the amplitudes N(ij,a) and N(ij,a) do 
not depend on signs e,e. 

We shall now show that for each set of the moduli i,j only one solution N(ij,a) and N(ij,a) 
is consistent with the remaining equations from the Table I. and that this solution resolves the 
four- fold sign ambiguity in relative phases of unnatural exchange amplitudes L and U. To this end 
we introduce a convenient notation for real and imaginary parts of reduced transversity amplitudes 
A = L,U,N 

A = A 1 +iA 2 = \A\ cos (xa + i\A\ sin a A (7) 
A = A\ + iA 2 = \A\ cos a>A + i\A\ sinciU 

where the phases oca = a^i, e) and Txa = aU(.7,e) for amplitudes A = L,U are given by (3.19). 
The equations for the remaining density matrix elements from the Table I. then take the form 

r 2 = >/2Im(p2)oiE = LiM + L 2 N 2 - L 1 N 1 - L 2 N 2 (8) 
s 2 = V2Im(p° z )l\-E = L 2 N ± - L ± N 2 + L 2 N 1 - L{N 2 

r 3 = -/m(^)_nS = FiJVi + U 2 N 2 - U 1 N 1 - U 2 N 2 (9) 
s 3 = -/m(^)-nS = C7 2 Aq - FiiV 2 + U 2 N 1 - U l N 2 

The equations (8) and (9) are four linear equations for four unknowns N\,N 2 , Ni,N 2 with four 
solutions Ni(ij,e,e), N 2 (ij, e,e), N\{ij,e,t) and N 2 (ij,e,e) corresponding to e = ±1 and e = ±1. 
One of these four solutions of (8) and (9) must be equal to one of the solutions of (1). To be specific, 
let us suppose that this physical solution for N(ij, a),N(ij, a) is the solution with the positive sign 
a = +1 in (4). It corresponds to some specific signs e + and e + so let us label all amplitudes in this 
solution A + ,A + . We now show that the amplitudes N~,N for the other solution with a = — 1 
cannot be a solution of (8) and (9) for any choice of signs e and e of the phases in (3.19). 

To prove this statetement let us suppose the contrary and assume that the amplitudes N~ , N 
satisfy (8) and (9) for some amplitudes A~ ,A , A = L,U from (3.19) corresponding to particular 
values of signs e~ and e~. Since the amplitudes A~,A , A = L,U differ from the amplitudes 
A + ,A + , A = L,U only in the sign of phases, their real parts are the same and the imaginary parts 
differ at most by sign 

A 2 =\A+, A; = \A+, A = L,U (10) 

where A = ±1,A = ±1. In the next step we subtract the set of equations with amplitudes N~,N 
from the set with amplitudes N + , N + and get a homogeneous set of equations 

Ll(N+ - N~) + Lt{N+ -AAV) - L+(Nt -NT) ~ L+(N$ - XN 2 ) = (11) 
L 2 (A^ - AAf ) - L^(N 2 + - N 2 ) + L 2 {Ni - XN j~) - Li(N 2 - N 2 ) = 

Ut(N+ - Af ) + Vt{N+ - \N 2 ) - U+(N+ - N~) - U+(N+ - XN 2 ) = (12) 
U 2 h (N+ - AAf) - F+(A+ - A 2 -) + U+(N+ - XN~) - *7+(A 2 + - A 2 ) = 
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We now write (6) in the form 

coso^ = -^j- (am T CL2S1) (13) 

, 1 , V 

sma N = (aisi ± a 2 r 1 ) 

and calculate differences and sums 

N+ - JVf = -2o2ai, JV+ + iVf = +2ain (14) 

N+ - N 2 = +2a 2 n, N+ + JVf = +2oiSi 

We can write similar expressions for amplitudes iV + and iV by replacing the paramerters a\ , a 2 in 
(14) with parameters a±, a 2 corresponding to the form (13) of the equations (5) and (4) for cos ajy 
and sin oat. 

We now examine the equations (11) and (12) for each possible choice of A and A. The case 
A = A = +1 is excluded as the amplitudes N~, N cannot satisfy the same system as the amplitudes 
N + , N + . For the case with A = A = — 1 the equations (11) and (12) take the form 

si(-L + a 2 + L^ai + L + a 2 - L^a^j = (15) 
rA -h\a 2 + -L + ax - Lfa 2 + L + a"i J = 



s^-Ut^ + ut^ + U+^-U+a^ = (16) 
n (-Utcb + Utax - U^a 2 + f/ 2 +oi) = 

The terms in the parentheses are all different. In particular, the large differences between the 
moduli \L\, \L\ and \U\, \ U\ mean large differences in the parentheses for r\ at least one of which 
must be non-zero. Since r\ has been measured in the CERN experiments on polarized targets and 
it is non-zero, the equations (15) and (16) cannot be satisfied and this case is excluded. 
For the case A = +1 and A = — 1 the equations (11) and (12) read 

{-h\a 2 + L 2 ai + Lfa^sx - ^L + a 2 ^rx = (17) 

(-L + a 2 +L 2 °i ~ L^a^rx - ^L + a 2 ^si = 

(-17+02+17+01 + f/ + o 2 )si - (ll+a^n = (18) 

(-f7 + a 2 + U 2 ai - £/ + a 2 )ri - (u 2 + a 2 js 1 = 

Combining the first two and the last two equations we obtain two equations of interest 

(-L + a 2 + L + ax + L + a 2 )si = (-L + a 2 + L + oi - L^a^rf (19) 
(-TJ + a2+TJ + ai + [/ + o 2 )s? = (-TJ 1 a 2 +I7 + oi-[/ 1 + 02)r 1 2 

The terms in the parentheses are all different. The parentheses for amplitudes L, L and amplitudes 
U, U are not proportional to each other since the phases of L, L and U, U are 180° out of phase, 
as seen in Fig. 7. Moreover, as we have discussed in Section VII., the moduli of these amplitudes 
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have different p°(770) widths and structures around /o(980). The equations (17) and (18) thus 
cannot be satisfied and this case is excluded. The analysis of the case A = — 1 and A = +1 is similar 
with the same conclusion. The solution N + , N + selects a unique set of phases of amplitudes L, U 
and L, U corresponding to e + and e + . The equations (8) and (9) change when the phases of these 
amplitudes change sign. The change results in a different solution for amplitudes N, N which is 
not compatible with the data in (1). 

We conclude that the measurements of S- and P-wave density matrix elements Imp® and Imp® 
unambiguously select a unique solution for S- and P-wave reduced transversity amplitudes. For any 
set i,j of solutions for the moduli there exists only one solution for the phases of natural exchange 
amplitudes and only one set of signs of phases of unnatural exchange amplitudes in (3.19). The 
corresponding solution of equations (8) and (9) for Ni,N 2 , N\, N 2 is the only solution that satisfies 
the conditions \N(i)\ 2 = N? + TVf and |A^(i)| 2 = + n\. 
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